Presence of energy flux in quantum spin chains: 
An experimental signature 
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Using the XXZ model for the description of one-dimensional magnetic materials we show that 
an energy flux, Jb, produces a shift, 8k ~ ■s/Je, in the characteristic wavenumber of the spin-spin 
correlations. We estimate 8k for a realistic experimental setup and find that it is measurable in 
inelastic neutron scattering experiments. 
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I. INTRODUCTION 

The problem of thermal transport in one-dimensional 
systems has been much investigated, the main goal being 
to derive Fourier's heat law. Analytical and numerical 
studies of a number of classical lattice-dynamical models 
indicate that the condition for Fourier's law to hold is 
the presence of strong nonlinearities i.e. nonintegr abil- 
ity (chaoticity) of the dynamics pj . Although quantum 
systems have been less studied, it appears that similar 
considerations apply to quantum spin chains as well j2j. 

Integrable systems, on the other hand, show anoma- 
lous thermal transport. No internal thermal gradient is 
formed in a harmonic crystal H or in a transverse Ising 
chain Q and, as a consequence, the energy (heat) flux is 
not proportional to the temperature gradient inside the 
sample. The origin of this anomaly may be the fact that 
the energy current in integrable systems often emerges as 
an integral of motion which automatically yields anoma- 
lous thermal transport coefficients H . 

The flat temperature profile in the presence of energy 
current is an intriguing feature of integrable systems, fn 
effect, it points to the existence of a homogeneous state 
carrying finite energy current. In this paper, we shall ex- 
plore the experimentally measurable properties of such a 
state by studying the XXZ spin chain in the presence of 
an energy current. The most spectacular feature of such a 
state is the incommensurability of magnetic excitations. 
Namely, in the presense of energy flow je, the charac- 
teristic wave vector is shifted from its antiferromagnetic 
value 7r by the amount 8k ~ \[Je- Since there are quite a 
few well established realizations of quasi-one-dimensional 
Heisenberg chains [§ (e.g. KCuF 3 §, Cs 2 CoCl 4 §, Cop- 
per Benzoate 0], Sr2Cu03, CS2CUCI4 [||), we believe 
that the predicted changes in the dynamical correlation 
functions bear direct experimental relevance. 

The basic problem of constructing a state which car- 
ries an energy current is the noncquilibrium nature of 
that state. Even if we assume that the flat temperature 
profile means the existence of equilibrium, we still face 
a problem that the value of the established temperature 
is not known M). We shall avoid this problem by re- 
stricting our calculation to zero temperature (T = 0) 



and assuming that the ground state correlations are ro- 
bust enough to survive at low temperatures. 

The construction of a homogeneous state with energy 
current at T = can be done by adding the energy cur- 
rent with a Lagrange multiplier to the XXZ Hamiltonian 
and then finding the ground state. Similar calculations 
have been carried out already for the transverse Ising 
and XX chains Ojll and, in a different context, for 
the XXZ model |fj,|lj. The new result we report is 
the calculation of an experimentally accessible parame- 
ter, namely the shift, 8k, of the characteristic wavenum- 
ber in the spin-spin correlations as a function of the en- 
ergy current, je- 

Once we have Sk(jE), we turn to a realistic experimen- 
tal setup and estimate je flowing through a single spin 
chain which gives an estimate of 8k. Our result shows 
that 8k is in the accessible range of an inelastic neutron 
scattering experiment. 



II. THE MODEL AND THE CHARACTERISTIC 
WAVE NUMBER 

The model we study is the spin- 1/2 XXZ chain defined 
by the Hamiltonian 



(1) 



where the spins af (a = x, y, z) are Pauli spin matri- 
ces at sites I = 1, 2, N of a one-dimensional periodic 
chain (cr% +1 = sf). We shall use the parametrization 
A = cos 7 and consider only the 'antiferromagnetic' re- 
gion 0<7<7r/2. In order to impose a fixed energy 
current je in the ground state, we add the current oper- 
ator to the Hamiltonian with a Lagrange multiplier 



H = iT x 



Ail 
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where 
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Importantly, ]e is an integral of motion, \Je, H\ = 0, 
thus indicating that i) the transport of energy is singular 
in this system Q and ii) the states carrying fixed energy 
current can be obtained as stationary states of H xxz . 

The XX Z model can be described in terms of interact- 
ing fermions and has been solved using the Bethe Ansatz 
method. The same approach works in the presence of 
the driving term, \]e, as well, and the solution has been 
given in 13 1%. An interesting feature of the solution 
is that the system displays rigidity against the drive, 
namely the ground state supports a nonzero energy cur- 
rent, {]e) = 3e 7^ 0, only if the coupling A exceeds some 
critical value A c (7). As we are interested in fixed energy 
currents, we simply choose sufficiently large values of A. 
Furthermore, since in realistic situations je turns out to 
be small, we concentrate on the region A « A c (7), in 
which case ]e oc (A — A c ). 

Once the energy current flows, an important restruc- 
turing takes place in the ground state. The single Fermi 
sea characterizing the ground state without current splits 
into two Fermi seas as shown in Fig.l for the simple case 
of the XX limit (A = 0) where a free-fermion description 
applies. There are now four Fermi wave vectors, ±7r/2 
and 7r/2 ± 8k, and the structure of the ground state im- 
mediately implies that there will be gapless excitations 
at wave vectors 0, 8k, 25k, tt — 8k, tx and n + 8k, which is 
readily confirmed by the exact solution at arbitrary A. 
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FIG. 1. Single-particle fermionic spectrum in the XX limit 
(A = 0) of the -Hxxz + Aj'b hamiltonian with and without en- 
ergy current in the ground state (dashed line, A = 1.3 and 
solid line, A = 0, respectively). The energy is measured in 
units of J while the wave number, k, is given in units of the 
inverse lattice spacing. The Fermi energy is zero indepen- 
dently of A. 



Thus an incommensurability characterized by 8k ap- 
pears in the system. This can be seen readily in the 
ground-state correlations. Indeed, it has been shown 
Jl2| , fi"5| that, in the XX limit, the longitudinal correla- 
tions for small je ^ can be expressed in a scaling form 



T e a e+n)jE=Q 



<&(8kn) 



(4) 



where &(x — * 0) = 1 and the large argument asymptotics 
of the scaling function is given by 



lim 



1 



(1 



(5) 



Since (of cr| +n )j B=0 — (-l) n /y/n equations (|,||) imply 
that, as the current is switched on, the static structure 
factor develops additional peaks at k = tt ± 8k (as it will 
turn out, 8k is small thus it is better to speak about the 
k = tt peak developing shoulders for j E ^ 0). 

In order to connect 8k to the current one determines 
both je and 8k through A and then eliminates the La- 
grange multiplier. The expressions arc simple for the XX 
limit [HH 



JE 



2tt% \ A 2 



cos 8k = X 



and, for small currents (A > A c = 1), they yield 



8k = 




(6) 



(7) 



3e^ 



= J 2 /h, has 



where a 'natural unit' of the current, 
been introduced. 

The above calculation can be carried out for any < 
A < 1 and the result for 8k differs only in a prefactor of 
order unity |l6[ 



8k 
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JE 



ir sin j V jW 



As one can see, the largest 8k is obtained in the XX limit 
(7-7T/2). 

In principle, if 8k is large enough then the extra peaks 
at 7r± 8k should be observable as Bragg peaks in an elas- 
tic neutron scattering experiment. In practice, however, 
the incommensurate modulations of distinct spin chains 
are not correlated and, as a consequence, the delta func- 
tion of the Bragg peak would spread out into a plain and 
the effect would be unobservable. 



III. STRUCTURE FACTOR AND AN 
EXPERIMENTAL SETUP 

It is more promising to look for an experimental signa- 
ture in an inelastic neutron scattering experiment where 
the excitations of the system are measured and no coher- 
ence among the chains is needed. Taking into account 
the facts that, for je = 0, most of the spectral weight is 
concentrated on the region around the antiferromagnetic 
wave vector 7r, and furthermore that, for ]e ^ 0, there 



2 



are gapless excitations at wave vectors ir±5k, we expect 
that the presence of the current manisfests itself via the 
emergence of additional inelastic peaks at wave vector 
7r±<5fc. This expectation can be put on a more solid base 
by calculating the dynamic structure factor and examin- 
ing the relative weights at wave vectors 7r and tt ± Sk. 

The simplest case is again the XX limit where the 
calculation of the time-dependent transverse correlation 
function, (o-*(t)<7o(0)), is straightforward. There is, how- 
ever, a principal problem at the outset of the calcula- 
tion. Namely, it is not clear whether the time-evolution 
of (Tn(t) is governed by H xx or by ff xx — Xj E . We shall 
take the view that in reality the current-carrying state is 
formed as a result of boundary conditions. Thus the local 
perturbation caused by an incoming neutron evolves by 
the local hamiltonian i.e. by i? xx [jiTj. Once (ct^(^)ctq(O)) 
is known its Fourier transform in time and space gives the 
structure factor S zz {k,uj) as shown in Fig. 2. 
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signature for the current-carrying state. The remaining 
question now is how to estimate 5k. 

As we can see from (||) , an estimate of Sk requires the 
value of the energy current, je- Thus we should, in prin- 
ciple, calculate je in a spin chain where the two ends 
are kept at different temperatures. We are unable to do 
this for any reasonable size system, and so we shall treat 
the energy flux as a parameter taken from experiments 
(Je = j E )• Then a thermodynamic measurement of 
j E xp can be used to estimate the value of Sk in an inde- 
pendent neutron-scattering experiment. Below we shall 
show how to estimate using parameters from a real- 
istic experimental setup. 
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FIG. 3. Experimental setup for measuring the effect of en- 
ergy flux. The solid lines represent the spin chains in a cubic 
sample of volume Z 3 with d being the distance between the 
chains. The energy flux is generated by keeping the two ends 
of the chains at temperatures T and T + 5T, respectively. 
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FIG. 2. The structure factors S zz (k,uj) displayed on the 
wavenumber-frequency (k — uj) plane for cases of (a) A = 1.00 
(no flux) and (b) A = 1.01. The unit of uj is J '/h while k is 
measured in units of the inverse lattice spacing. The darkness 
of the shading is proportional to S zz (k, uj). 

As we can see, a large part of the weight of the jE = 
peak of the structure factor at ir shifts to tt ± 5k for 
Je 7^ 0. Thus one can expect that even if 5k is small, the 
presence of a small energy current will result in a broaden- 
ing by 25k of the inelastic peak centered at wavevector tt. 
It is this broadening that we propose as an experimental 



Let the sample be a cube of side I = 10 _2 m and let 
the spin chains be along x direction with the distance be- 
tween the neighboring chains being d = 10~ 9 m. Further- 
more, let the sides of the cube perpendicular to the chains 
be at temperatures T and T+5T (see Fig. 3). The temper- 
ature should be chosen to be low in order to minimize the 
phonon contribution to ■ However, T cannot be too 
small since then the small coupling between the chains 
makes the system three-dimensional. The value of 5T 
should be, in principle, chosen large but the limitations 
of cryogenics of a realistic setup restrict the steady-state 
temperature-differences to ST < OAT. From the above 
considerations we arrive at the following ranges for the 
possible temperatures and temperature-differences 

T = (1-10)°K ; ST = OAT = (0.1 - 1)°K . (9) 

The total current of heat across the sample can now be 
estimated as 



■total 
JE 



(10) 



provided we know the heat conductivity, k. We note here 
that the finiteness of the experimental k is not in contra- 
diction with the singular nature (n mt = oo) of the heat 
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conductivity of integrable spin chains. A macroscopic 
sample consist of spin chains of characteristic length 
I « 10~ 4 cm and the energy must also be transported 
between chains. This leads to the loss of ideal conductiv- 
ity and results in a finite k. Consequently, the estimate 
of energy flux using the experimental k does give an esti- 
mate of the energy flux through the chains provided the 
spin chains are the main channels of energy transport. 

Unfortunately, k, is not available for the materials we 
have in mind, |^-|^], and another problematic issue is how 
much of the conductivity comes from the spin-chains. 
Since measurements of the magnetothermal conductiv- 
ity of magnetic materials Jl8|,|l9| indicate that spin waves 
provide a significant fraction of the low-temperature ther- 
mal conductivity, we shall assume that an order of mag- 
nitude estimate of the energy current through the spin 
chains is given by jg . Furthermore, we shall assume 
that, as a value of k, we can take a characteristic value 
of this parameter in crystalline magnetic materials in the 
temperature range T = (1 - 10)°^ f§: 
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it«(l- 10) 



w 



m ■ °K 



(11) 



We can then estimate j^ tal w (10~ 3 - VQ- l )W and, 
since the number of spin-chains in the sample is J\f = 
l 2 /d 2 = (10- 2 /10~ 9 ) 2 = 10 14 , we obtain the energy flux 

-\ ■ -exp 

per chain, j E , as 



■ exp 

Je 



■total 

N 



= (10~ 17 - \Q- lb )W 



15 \ 



(12) 



As we have seen (^) the natural unit of energy current 
in a spin chain is j$ = J 2 /h. Using a characteristic 
value of J = (1 — 10)°K for the spin coupling we find 
w (10~ 12 - W- 10 )W and obtain the following esti- 
mate for the shift of the wavenumber 




10" 



(13) 



This is our central result. Since Sk ~ 10~ 2 is accessible 
in an inelastic neutron scattering experiment, the effect 
of shift in the wavenumber should be observable. 

In summary, we have studied an integrable system 
which doesn't obey Fourier's law. We proposed that, un- 
der some simplifying assumptions, one can explore states 
of this system which carry current of energy and, fur- 
thermore, one can derive theoretical results verifiable in 
experiments. 
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